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C^ , Abstract 

The spectrum profile that emerges in molecular spectroscopy and 
atmospheric radiative transfer as the combined effect of Doppler and 
pressure broadenings is known as the Voigt profile function. Because of its 
convolution integral representation, the Voigt profile can be interpreted 
as the probability density function of the sum of two independent random 
variables with Gaussian density (due to the Doppler effect) and Lorentzian 
*vS ', density (due to the pressure effect). Since these densities belong to the 

^^ . class of symmetric Levy stable distributions, a probabilistic generalization 

is proposed as the convolution of two arbitrary symmetric Levy densities. 
We study the case when the widths of the considered distributions depend 
1^^^ ■ on a scale-factor r that is representative of spatial inhomogeneity or 

(«— «, ' temporal non-stationarity. The evolution equations for this probabilistic 

generalization of the Voigt function are here introduced and interpreted 
as generalized diffusion equations containing two Riesz space-fractional 
derivatives, thus classified as space-fractional diffusion equations of double 
order. 
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1 Introduction 

In molecular spectroscopy and atmospheric radiative transfer, the combined 
effects of Doppler and pressure broadenings lead to the Voigt profile function, 
which turns out to be the convolution of the Gaussian (due to the Doppler 
broadening) and the Lorentzian (due to the pressure broadening) distributions. 
The study of the Voigt profile is an old issue in literature and many efforts have 
been directed to analyze its mathematical properties and relations with other 
special functions and to obtain its numerical computation, e.g. [1, 2, 4, 5, 7, 9, 
10, 12, 16, 17]. 



In most papers appeared in the literature, the Voigt profile is defined in terms 
of a single weight-parameter, that is the ratio of the Lorentzian to Gaussian 
width. In physical applications, this parameter indicates which distribution is 
more important than the other. Generally, the weight-parameter is considered 
with a constant value fixed by the process. Here we are interested to study 
the Voigt profile function when the widths are not constant but depending on 
a scale-factor with a power law. Physically, the one-dimensional variable of the 
Voigt function is a wave-number and then this permits to take into account 
spatial inhomogeneity or temporal non-stationarity when the scale-factor is the 
distance from an origin or the elapsed time from an initial instant, respectively. 

Since in probability theory the Gaussian and the Lorentzian distributions 
are known to belong to the class of symmetric Levy stable distributions, 
in this framework we propose to generalize the Voigt function by adopting 
the convolution of two arbitrary symmetric Levy distributions. Moreover, 
we provide the intcgro-differential equations with respect to the scale-factor 
satisfied by the generalized Voigt profiles. These evolution equations can be 
interpreted as space-fractional diffusion equations of double order. 

The paper is organized as follows. In section 2 the basic definitions for the 
Voigt profile are given. In section 3 the connection with the class of Levy stable 
distributions and the probabilistic generalization are introduced. In section 4 
we derive the integro-differential evolution equations of the generalized Voigt 
function with respect to the scale-factor. In section 5 the limits of low and 
high scale-factor values are investigated. Finally, section 6 is devoted to the 
concluding remarks. 

2 The Voigt profile function 

The Gaussian G{x) and the Lorentzian N{x) profiles are defined as 



G(x) = —= — exp 
Vttwg 
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where log and lol are the corresponding widths. From their convolution we have 
the ordinary Voigt profile V{x) 

r+°° LOT /ujr /■+°° c^ii/'^Gf 



The main parameter of the Voigt function is the weight-parameter a defined as 
a = lol/loq^ which is the ratio of the Lorentzian to Gaussian width and then a 
measure of the relative importance between their infiuence on the properties of 
the process. Generally, the a < 1 case is important in astrophysics while a > 1 
in spectroscopy of cold and dense plasma [4]. In particular, two limits can be 
considered: i) a ^ 0; ii) a — > oo. In the first case, the Lorentzian contribution is 
negligible in respect of the Gaussian one, while the inverse occurs in the second 
case, see Fig. 1. 



a=0.01 



a=0.1 



0.001 




0.01 



1 



0.1 



0.01 



0.001 



^N. 


1 




1 

Voigt 

Gauss 


1 

------- 


- % 1 






Lorentz 


H 


: Q 


^ A. 






: 






- 


1 




k 




>^ 




0.01 



1 



Figure 1: Comparison between the Voigt, Gauss and Lorentz distributions with 
weight-parameter a ~ lol/loq — 0.01, 0.1, 1, 2. 



Let /(k) be the characteristic function, which is the Fourier transform of /(a;), 
so that 
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then V{k) = G{k)N{k) — e '^g'* /^e '^^1'*!. These formulae imply the following 
integral representation for the Voigt profile 
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V{xJ^—l e-*'*^e-"^'''/''-"^l'=ldK= - / c-'^^'^-'^'^^'/^ cos{Kx)dK. (4) 



3 The probabilistic generalization of the Voigt 
profile function 

It is well known that if Xi and X2 are two independent random variables with 
probability density function (PDF) qi and q2, respectively, then the PDF p{w) 
of the random variable W = Xi + X2 is given by the convolution integral 



p{w) 
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qi{w - X2)q2{x2) dx2 



(5) 



From (2) and (5), the Voigt profile can be seen as the resulting PDF of the sum 
of two independent random variables, one with Gaussian PDF and the other 
with Lorentzian PDF. 



The Voigt function has been generahzed in hterature in different ways, e.g. 
in [12] the integrand function in formula (4) is multiplied for a polynomial, 
in [14, 15, 6] the cosine function in (4) is replaced with the Besscl function 
and with the Wright function with one and more variables. In this paper we 
propose a probabilistic generalization in the framework of Levy distributions. It 
is well known that the Gaussian and the Lorentzian distributions are two special 
cases of the class {La{x)} of the symmetric Levy stable distributions, where a, 
< a < 2, is called characteristic exponent. Writing the characteristic function 
(3) of Levy distributions as Lain) — e~''^' , the Gaussian and the Lorentzian 
profiles defined in (1) are recovered with a = 2 and luq — 2 and with a = I and 
(jJl = ^, respectively, 

G{x)^L2{x)^-^e-^'/\ 7V(x)=ii(x) = i^i-. (6) 
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The Voigt function can be straightforwardly generalized in probabilistic sense by 
considering the sum of two independent random variables with symmetric stable 
densities. Mathematically, this corresponds to the convolution of two arbitrary 
symmetric Levy densities of characteristic exponents ai and a2- Denoting 
with V(x) the generalized Voigt function, its integral representation and its 
characteristic function V{n) are 

V{x)^ L^,{x~i)L^Md^, V(«) = e-I«l''^-I'^l°\ (7) 



4 The evolution equations with respect to the 
scale-factor 

In the previous sections the widths ujg and lol are considered to be constants 
and the weight-parameter a fixed. However, differently from the previous 
papers on the topic, we would like to know what happens when the widths 
loq and lol change in space or time with a power law with respect to a scale- 
factor. This analysis is relevant in inhomogeneous or not stationary cases, for 
which the scale-factor corresponds to the distance from an origin or the elapsed 
time, respectively. Conversely, constant values of widths can be considered for 
homogeneous and stationary cases. In the present section we consider a scale- 
factor T for both spatial inhomogencity and temporal non-stationarity. 

It is well known that the Levy density functions La{x,T) are the fundamental 
solutions of the space-fractional diffusion equation [3, 8, 13] 

^h^pll^D^L^ix,T), L„(x,0) = <5(x), 0<a<2, (8) 

where D" is a pseudo-differential operator known as the Riesz space-fractional 
derivative of order a. Such pseudo-differential operator is defined in terms of 
its symbol — |k|", i.e. the Fourier transform of D"f{x) is — |k|"/(k). We recall 



the explicit representations for a ^ 2 

( r(l + a)^. (a^^ f+^ fjx + 0- 2,f{x) + f{x - Q ^^ 
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(9) 
and the hmit D"f{x) ~ (Pf/dx^ when a ~ 2. The solutions of (8) have the 
power law scaling property 

L„(x,r)=r-V"L„(-^) . (10) 

For analytical and graphical representations of stable densities we refer the 
reader to [8]. We observe that the distributions defined in (10) are self-similar 
and they obey to the same power law scaling for all values of the scale-factor t. 
In particular, for a = 2 and a = 1 the Gaussian and the Lorentzian densities 
are recovered, respectively, and from (10) we have loq cx t^' ^ and w^ ex t. 

Following scaling (10), our generalized Voigt function (7) becomes 

V(x,r)=r-V"i-i/". r°^L^^ f^V^. { Jr-] di , (11) 
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and its characteristic function is 

V(k,t) ^e^l'^l"'^-!'^!'"^. (12) 

In this case, it is possible to show that formula (11) is the solution of the 
following integro-differential equation 

dV 

— ^ Dt'V{x,T) + DZ'V(x,t) , V(x,0)^5{x). (13) 

In fact, after the Fourier transformation Eq. (13) becomes 

|^-~|Kr^V(«,r)-|Kr^V(«:,T), V(«,0)-1, (14) 

which is solved by (12). When ai — \ and a2 ~ 2 the integro-differential 
equation (13) is the evolution equation of the ordinary Voigt function. Equation 
(13) can be classified as space- fractional diffusion equation of double order. The 
evolution of V(x, t) for different pairs of (ai, a2) with t = 0.1, 1, 10 is shown in 
Fig. 2. 

5 The asymptotic scaling laws for low and high 
scale-factor 

The Voigt (2) and the generalized Voigt (7) profiles are derived from the 
convolutions of two self-similar processes with different scaling laws and, as 
a consequence, the similarity is lost. However, we ask which are the scaling laws 
of the Voigt functions in the limits of low and high values of the scale-factor r. 
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Figure 2: The evolution of the generahzed Voigt function V(x,t) for the pairs 
(ai, aa) = {(0.5, 2), (1.5, 2), (0.5, 1.5), (1, 1.5)} with t ^ 0.1, 1, 10. 



Since for Levy stable densities with a ^ 2 the mean square displacement 
diverges, the same occurs for the ordinary and the generalized Voigt functions. 
Then, to analyze the scaling laws when r ^ and t —^ oo the variance 
{x"^) cannot be used. However, the scaling law of a Levy stable process with 
characteristic exponent a can be studied by mean the quantity {\x\'^)^''^, with 
< q < a. In this respect, we have 

(|a;|«) = 2/ x«V(a;,T)dx = 2V*(<7 + l,T), < g < min{ai, as} , (15) 
where V*(s) is the Mellin transform of V(x), x > 0, defined as [11] 



V*(s)=/ V{x)x'-'^dx, V{x) = —- V*(s)a 



' ds . 
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Without lost in generality, let us state ai < a2- Starting from (7) and following 
[10], the Mellin-Barnes integral representation of the generalized Voigt function 
V{x) is 



V(x,t) 
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Hence its Mellin transform is 

^(so-l)/a2 

V*(si,t) = r(so) cos(so7r/2) 
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and finally, 



(l^r) = 2V*(g+f,T 
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Applying the residue theorem to T I '* ^"^''^ ) ^ we obtain the convergent series 

for r ^ oo, 

(|.r) = -^^r(,+ f)sin(W2)f: ^^^ f^^^!^^) .-«(".M-i) . (20) 

n— ^ ^ 

Applying the residue theorem to r(si), we obtain the convergent series for 

T^ 0, 

(IxH = -^l!^r(, + 1) sin(W2) V ^^r (^ll^\ r"(--/-) . (21) 
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In papers [3, 13] the limits (r ^ 0,t ^ oo) are computed using a different 
method but with the same results. Then the two limits under consideration 



give 

(a;9)i/9 ocT^/"^ , T^O 
r ^ oo 
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In the limit r ^ 0, the corresponding scaling law of the generalized Voigt profile 
is governed by the Levy density with the higher value of the characteristic 
exponent, while in the limit r — > oo by that with the lower value. In particular, 
for the ordinary Voigt profile (ai = 1 , q;2 = 2) the process scales as r^ " and r 
for law and high values of the scale-factor, respectively. This means that if the 
power law represents inhomogeneity, or non-stationarity, the resulting profile 
is approximated by a Gaussian for small distances from an origin, or small 
elapsed times, and it is approximated by a Lorentzian for large distances, or 
large elapsed times. This result is consistent with the usual limits a ^ and 
a -^ oo, see Fig. 1. 



6 Conclusions 

In the present paper we have considered the Voigt profile function and proposed 
a probabihstic generahzation as the convolution of two arbitrary symmetric 
Levy densities. Generally, the Voigt profile characteristics are studied with 
respect to a weight-parameter a that is the ratio of Lorentzian to Gaussian 
widths, a = uj[^/uja, and it is assumed to be a constant property of the process. 
Conversely, here we have considered both widths depending on a scale-factor r 
that is representative of inhomogeneity or non-stationarity. We have introduced 
parametric integro-difFerential equations for the ordinary and the generalized 
Voigt functions. These integro-differential equations can be classified as space- 
fractional diffusion equations of double order because they include two Riesz 
space-fractional derivative of different orders. In this respect, the present 
paper shows an application in physics of the distributed fractional derivatives 
formalism. 

Finally, the limits of the Voigt function for low and high values of the scale- 
factor are considered. The Voigt function turns out to be not self-similar, even 
if it is expressed as the convolution of two self-similar Levy processes. Its scaling 
law is dominated by the Levy density with the higher value of the characteristic 
exponent when t ^ and by that with the lower value when t ^ oo. 
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